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ABSTRACT 

We present new late-time near-infrared imaging of the site of the nearby core-co llapse super- 
nova SN 2012aw, confirming the disappearance of the point source identified by [Fraser et al. 
( |2012| ) and [Van Dyk et al.| ( |2012| ) as a candidate progenitor in both J and Ks filters. We re¬ 
measure the progenitor photometry, and find that both the J and Ks magnitudes of the source 
are consistent with those quoted in the literature. We also recover a marginal detection of the 
progenitor in F^-band, for which we measure F^=19.67 ± 0.40 mag. Comparing the luminosity 
of the progenitor to stellar evolutionary models, SN 2012aw appears to have resulted from the 
explosion of a 12.5±1.5 Mq red supergiant. 
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1 INTRODUCTION 


Core collapse supernovae (SNe) are luminous and energetic tran¬ 
sients associated with the demise of an evolved massive star. They 
play a critical role in the production and distribution of metals 
within galaxies (e.g. [Nomoto et al.|2006] l, and regulating star for¬ 
mation and galaxy evolution through feedback (e.g.|Scannapieco et| 
|al.|2008| ). Nearby supernovae with identified progenitors can also 
be used to test the predictions of stellar evolutionary models and un¬ 
derstand the connection between the final stages of massive stellar 
evolution, and the resulting explosion and compact remnant (|Heger| 
let al.|2003[|Eldridge et al.|2008] l. 

In the last decade the detection of progenitor candidates for 
nearby core-collapse supernovae in archival imaging has become 
relatively routine (see |Smartt||2015| for a recent review). To con¬ 
firm these candidates requires deep images to be taken after the SN 
has faded, in which the progenitor candidate should be no longer 
visible. This technique has been used to verify the progenitors of 
some nearby SNe < |Van Dyk|2013[|Maund et al.|2014||, but has also 
revealed several erroneous identifications ( |Maund et al.|2015| ). 

SN 2012aw is a nearby, well studied Type Il-Plateau (IIP) SN 
in M95, discovered by [Fagotti et al.| ( [20T^ . Followup observations 
of SN 2012aw were presented by [Bose et al.[([2013|) and [Dali’Ora[ 
leTaLlpOM] ); the SN appears relatively bright and quite si milar to 
SN 1999em, with an ejected ^®Ni mass of 0.06 M©. Both [Fraser 
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let al.| ( [20T^ and I Van Dyk et al.H2012^ (henceforth |F12| and 
identified a progenitor candidate for SN 2012aw in archival Hub¬ 
ble Space Telescope and ground based near-infrared (NIR) data. 
The candidate had an apparent magnitude in F814W (~/-band) of 
23.30±0.02 mag, however its colour was redder than one might 
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expect for even the coolest red supergiants, with F555W—F814W 
between 3.1 and 3.3 mag, leading both [F12[ and [vD12| to suggest 
that the progenitor suffered from significant extinction. [vD12[ and 
[F12[ estimated that the progenitor had a zero-age main sequence 
(ZAMS) mass of 17-18 M© or 14-26 M© respectively. However, 
[Kochanek et al.[ ( [2012[ ) pointed out that both of these mass ranges 
were likely overestimates due to an incorrect treatment of the ef¬ 
fects of circumstellar dust. [Kochanek et al.[ argued that a purely 
silicate dust composition would be more appropriate for a red su¬ 
pergiant, and that dust emission in the NIR must be considered. 
Moreover, they point out that scattered photons from a shell of dust 
around a star will also contribute to its observed flux; resulting in 
a tendency to over-estimate the effects of extinction on luminosity 
when using a [Cardelli et al.[p^8^ extinction law. 

In this Letter, we present new NIR observations obtained three 
years after the explosion of SN 2012aw, which we use to verify the 
identiflcation of the progenitor candidate found by both [F12[ and 
[vD12[ We adopt a distance of 9.8±0.2 Mpc to M95 (based on the 
average of the [Freedman et al.[2001[ and [Rizzi et al.[2007[ distances, 
derived from Cepheids and the Tip of the Red Giant Branch respec¬ 
tively; these are consistent with the 10 Mpc distance used by [F12| l. 
We take the foreground (Milky Way) extinction to be Ay =0.08 
mag; consistent with values used in [F12| while in in Sect. [^ we 
adopt the values for circumstellar extinction around the progenitor 
from [Kochanek et al.[ ^20l2 \ . 


2 OBSERVATIONS AND DATA REDUCTION 

New observations of the site of SN 2012aw were obtained with the 
8.2 m Very Large Telescope (VLT) + High Acuity Wide held K- 
band Imager (HAWK-I; [Kissler-Patig et al.|2008j l on the nights of 
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2015 March 2 and 7. HAWK-I consists of four 2048x2048 pixel 
detectors, with a combined field of view of 7.5'x7.5', and a pixel 
scale of 0.106"/pix. On March 2, M95 was observed using the Ks 
filter on MJD 57084.10, at an airmass between 2.0 and 1.6. On 
March 7, observations were taken on 57089.2 using the J filter at 
an airmass of between 1.2 and 1.3. 

For both J and Ks, four consecutive frames were taken at a 
single pointing, with the site of SN 2012aw approximately centred 
on the Q1 chip of HAWK-I, before the telescope was dithered ~5" 
to a new position. The site of SN 2012aw is ~2' from the centre 
of M95, and the background at this location is sufficiently sparse 
that separate off-target sky frames are not required. Each frame 
consisted of 3xl0s integrations (DIT=10 s, NDIT=3 in ESO ter¬ 
minology), which were read out in non-destructive readout mode. 
A total of 60 frames were obtained in J and 68 frames in Ks, to give 
combined on-source exposure times of 1800 s and 2040 s in / and 
Ks respectively. 

All data were reduced using the HAWK-I pipeline, running 
under the ESOREX environment. A bad pixel mask was constructed 
by identifying outlier pixels in dark frames. Twilight fiat fields for 
J and Ks were created, and were used along with the bad pixel 
masks to reduce the science frames. The most critical step in the 
reduction of NIR images is the removal of the sky background; for 
our data we employed a two-step procedure. A first pass sky im¬ 
age for each frame was created by median combining the 8 images 
taken at the two preceding dither positions, and the 8 images taken 
at the two subsequent dither positions. The sky image associated 
with each frame was then subtracted to give the first pass reduced 
images, which were then shifted and co-added. The combined first 
pass image was then used to identify point and extended sources, 
and create a mask for these. This mask was then shifted as neces¬ 
sary and applied to each frame. The masked images were then used 
to create an improved sky frame for each image (again using 16 im¬ 
ages taken before and after each frame at different dither positions). 
The refined sky frame was subtracted from each of the original fiat- 
fielded, bad pixel-masked images to give the second pass reduced 
images. Einally, the second pass reduced images were corrected for 
geometric distortion using tabulated distortion coefficients, before 
being combined to give the final reduced image. 

The measured seeing of both the reduced J and Ks images 
was 4.6 pixels (0.5"). According to the ESO Ambient Conditions 
Database, both nights were photometric at the time the observations 
were taken. 

The pre-explosion images used in |F12| were re-reduced for this 
work. We considered both Js-band images from VLT -i- ISAAC 
taken on 2000 March 25, and the New Technology Telescope -i- 
SOFI H and Ks images taken on 2002 March 24. We note that the 
//-band images were not analysed in |F12| as the exposure times are 
relatively short, however, a careful reduction reveals a source at the 
location of SN 2012aw. The ISAAC data have a pixel scale of 0.15 
pix arcsec”^, while SOFI has a scale of 0.29 pix arcsec”^. 

For the SOFI data, a bad pixel mask appropriate to the epoch 
of the observations was obtained from the SOFI webpages. The 
reductions were carried out with IRAE, with the exception of de¬ 
termining the illumination correction frame, for which an ESOREX 
script was used. Flat field images for both H and Ks were derived 
using a set of dome flats taken both with and without lamp illumi¬ 
nation. As the dome fiat does not exactly match the illumination of 
the sky, an illumination correction frame was derived from obser¬ 
vations of a bright star which was imaged at 16 different positions 
on the detector. The bad pixel mask, flat field and illumination cor¬ 
rection were then applied to all science frames. A sky frame for 



Figure 1. Filter transmission curves for ISAAC Js and HAWK-I /; SOFI //; 
and SOFI Ks and HAWK-I Ks 

each filter was created from the mean of the science frames, us¬ 
ing sigma clipping to remove sources. As the images are relatively 
shallow, it was not found to be necessary to use a two pass sky 
subtraction. Finally, the images taken in each filter were shifted to 
a common pixel grid and co-added. 10 K-band images were used 
(each of which comprised 10x6 s DIT) for a total exposure time of 
600 s, while in //-band 4 frames were taken with 3 x 15 s DIT) for 
180 s on source. 

The ISAAC A-band data were reduced within ESOREX. Flat 
fields were derived from twilight sky observations, while a bad 
pixel mask was constructed from dark frames; these were then ap¬ 
plied to all science frames. Separate off-target sky frames were 
available for these observations, and so we used the off-source 
frames taken immediately prior and subsequent to each science 
frame to create a sky frame, which was then subtracted from the 
on-source images. The reduced, sky-subtracted images were then 
aligned within IRAE and co-added to produce a single, deep image. 
The location of SN 2012aw lies very close to the edge of the de¬ 
tector in some of the frames, and so we did not include these in the 
final stacked image. The combined image is comprised of 5 expo¬ 
sures, each of which has 4x30 s DIT, to give a total time exposure 
of 600 s. 


3 ANALYSIS 

Ideally, the late-time observations of the site of SN 2012aw would 
be taken with the same instrument and filter as the pre-explosion 
data. As this was not possible, we tested the effects of mis-matched 
filters using synthetic photometry. We obtained the throughput 
curves for the filters used from the ESO webpages, these are shown 
in Eig. We used the throughput curves together with |Castelli| 
|& Kuru^ ( |2004| ) models appropriate to a massive red supergiant 
and the PYSYNPHOT package to calculate the expected Jsisaac — 
7hawk-i and K^sofi—K^hawk-i zeropoint differences (effec¬ 
tively a colour) for the progenitor. Eor J, the largest difference is 
0.09 mag, while for Ks the difference is < —0.01 mag. 

To compare the pre-explosion and late-time images, it is nec¬ 
essary to register them to a common pixel grid. The pixel coordi¬ 
nates of point sources common to both frames were measured, and 
IRAE GEOMAP was then used to derive a geometric transformation 
between the matched coordinate lists. This transformation was then 
applied to the late-time image to map it onto the same pixel grid as 
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the pre-explosion frame. By transforming the late-time rather than 
the pre-explosion image, we avoid reducing the signal-to-noise of 
the progenitor candidate. We used the High Order Transform 
OF PSF And Template Subtraction (hotpants) package to 
match the Point Spread Functions (PSFs) of the two aligned frames 
and perform difference imaging between them. The resulting sub¬ 
tracted images are shown in Fig.[^ In both filters, the source visible 
in pre-explosion imaging has disappeared at late times, and appears 
as a source in the difference images, confirming the identification 
of the progenitor candidate. In the 7-band there is a residual cor¬ 
responding to a poorly subtracted nearby bright star, and there is 
some structure to the background in the subtracted image which is 
caused by proximity to the edge of the image, while in TT-band the 
subtraction appears good, with only a single source (most likely 
a variable foreground star) in the difference image apart from the 
progenitor. 

We used eight bright Two-Micron All Sky Survey (2MASS; 
[Skrutskie et al.|[2006] ) sources within the footprint of the entire 
HAWK-I mosaic to determine the photometric zeropoint for the 
HAWK-I Ks-hmid image, with an uncertainty in zeropoint of 0.11 
mag. These were in turn used to define a local tertiary photometric 
sequence of intermediate-magnitude sources which were visible in 
the pre-explosion SOFI Ks image. The same procedure was used to 
determine the zeropoint of the 7-band image, with an uncertainty of 
0.15 mag. We performed PSF-fitting photometry on the progenitor 
in template subtracted images in Ks and on the un-subtracted image 
in 7 using the SNOOPY packag^ and calibrated the resulting mag¬ 
nitudes using the local sequence. The magnitude of the progenitor 
was found to be 7=20.83 zb 0.18, and 7^5=19.56 zb 0.29 mag, where 
the error refiects both the uncertainty in zeropoint, and progenitor 
photometry as determined by artificial star tests. These magnitudes 
are consistent with the values found by both | vD 1 2| and |F 1 2] As an 
additional test, we also performed aperture photometry on the pro¬ 
genitor in the TT-band image (where the nearby source to the south 
east is not visible in the pre-explosion frame) and find a magnitude 
which is consistent to within the uncertainties with the results from 
PSF-fitting photometry. 

In principle, one could improve the photometry for the 7-band 
pre-explosion data by measuring the magnitude of the progenitor 
on the template subtracted image, which may reduce the contribu¬ 
tion of the nearby source to the south west to the measured fiux. 
However, as we were unable to obtain a satisfactory subtraction for 
the 7-band image due to the proximity of SN 2012aw to the edge of 
the detector in the pre-explosion image, measuring its magnitude 
on the difference image did not improve our results. 

As no late-time images were taken in the //-band, we were un¬ 
able to perform template subtraction, or to improve the photometric 
calibration of the pre-explosion data using observations of standard 
fields. We performed PSF-fitting photometry on the SOFI H image 
(shown in Fig. using SNOOPY, and estimated the photometric 
error through artificial star tests at the location of the progenitor. 
The zeropoint of the image was determined through aperture pho¬ 
tometry of six sources in the field for which 2MASS photometry 
were available. While 2MASS has a large photometric uncertainty 
for faint sources, the uncertainty in zeropoint is insignificant com¬ 
pared to the photometric error for the progenitor. We determine the 
magnitude of the progenitor to be 7/ = 19.67 zb 0.40. 


^ SNOOPY in an IRAF package based on daophot, and developed by. E. 
Cappellaro for SN photometry. 



Figure 3. Pre-explosion SOFI //-band image showing the detection of the 
progenitor. Inset in the upper right is a magnified view of the position of 
SN 2012aw covering 5" x5", with the progenitor circled. 


4 DISCUSSION 


Our photometry is consistent with that measured by|F12|and vD12| 
and so the progenitor luminosity calculated by |Kochanek et al.| 
( |2012| ) remains valid. We adopt the la range for the progenitor 
luminosity assuming silicate dust from |Kochanek et al.| who con¬ 
strain the progenitor to have a luminosity of log L/L0=4.8-5.O dex. 
Using the solar metallicity STARS models ( [Eldridge & Tout|2004] ) 
as presented in [Smartt et al.| ( |2009| ), this luminosity corresponds 
to a progenitor zero-age main sequence (ZAMS) mass of between 
11 and 14 M©. We hence adopt 12.5zbl.5 M© for the mass of the 
progenitor of SN 2012aw. We note that other modern stellar evo¬ 
lutionary codes will give similar results, as the final model lumi¬ 
nosities for a given mass are comparable to those from the STARS 
code( |Smartt et al.|2009| ) . 

I Van Dyk et al.| ( [20T5] ) presented late-time ultraviolet and opti¬ 
cal imaging from HST taken on 2014 April 23, in which a resolved 
light echo was observed for SN 2012aw, caused by diffuse inter¬ 
stellar dust scattering the SN flux. The SN itself was still visible in 
these observations as a point source, although rVan Dyk et al.| ( |2015) ) 
found that its magnitude in F814W was 0.4 mag fainter than that of 
the progenitor, and suggested that the latter had vanished. The light 
echo is more pronounced in the bluer bands, and hence will not af¬ 
fect our NIR observations (which were also taken ~1 yr after the 
7/57 data). 

The disappearance of the progenitor of SN 2012aw brings the 
sample of Type IIP SNe for which a red supergiant progenitor can¬ 
didate has been confirmed through late-time imaging to seven. SN 
2012aw is among the highest mass progenitors seen to disappear 
thus far, and despite some erroneous early claims that it may lie 
above the observed 16.5 M© limit for Type IIP SN progenitors, 
it now appears to be securely below this limit. In fact, even if 
we compare the la upper limit on the luminosity of SN 2012aw 
(log L/L©=5.0 dex), to the model luminosity at the end of core He- 
burning (which is a conservative lower limit on the luminosity at 
the point of core-collapse), then the progenitor must be <17 M©. 
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(d) SOFI Ks pre-explosion (e) HAWK-I Ks late-time (f) Subtracted image {Ks) 


Figure 2. The pre-explosion and late-time J and /sT-band images of the site of SN 2012aw, together with the subtracted images. Tick marks indicate the position 
of SN 2012aw, the scale and orientation of all panels are shown in the upper left. The inset in each panel shows a 5" x5"region centred on the SN location. In 
the inset, the late-time image has been transformed to match the pixel scale and orientation of the pre-explosion image. 


Two independent techniques have been applied to estimate 
the ZAMS mass of the progenitor of SN 2012aw. [Jerkstrand et al.| 
( |2014| ) used nLTE radiative transfer models to derive a progenitor 
mass from the observed strength of emission lines in the late-time 
nebular spectra of the SN. Based on this, |Jerkstrand et aT] ruled out a 
progenitor ZAMS mass >20 M© and found that a ~15 M© model 
consistently reproduced the observations. However, the 12 M© 
model of |Jerkstrand et al.| also provides a reasonable match to the 
line strengths, in particular for the sensitive [O l] AA 6300,6364 A 
lines. |Dair Ora et aL] ( |2014| ) used a radiation-hydrodynamics code 
to model the early luminosity, temperature and velocity evolution 
of SN 2012aw, and found a best fitting progenitor model had a 
19.6 M© envelope mass. This is clearly inconsistent with the pro¬ 
genitor mass from direct imaging, and moreover continues the ten¬ 
dency for hydrodynamic models to yield higher progenitor masses 
when compared to direct detections ( [Utrobin & Chugai|2008| l. We 
note that |DairOra et al.| { |2014| ) performed an initial exploration 
of the potential progenitor parameter space using a semi-analytic 
code, and that the distribution of their models as a function of 
mass appears to show a local minimum at ~16M©. However, the 
full radiation-hydrodynamics code was unable to simultaneously fit 
the plateau duration and velocities for a model with this mass. 
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